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ABSTRACT 

A method for predicting swell from tropical cyclones 
using a spectral wave model (TYWAVES) was tested. The model 
was applied to predicting swell propagating from three 
typhoons in the Western North Pacific through gaps in the 
Ryukyu Islands into a region of the East China Sea. The 
model involves a source region concept which considers 
only the swell emanating from regions of peak energy in 
moving typhoons. For three representative typhoons, 
predicted heights were not significantly different from 
the observed heights. The time of occurrence of the pre- 
dicted peak height agreed well with observational values 
for the swell from two typhoons, but lagged by 6-12 hours 
meme che third. 

The dominant swell period and direction predicted 
by the model were not verifiable by data available for 
this study. 

Shoaling and refraction effects were considered in the 
prediction, in a simplified way, but attenuation was 
ignored even for the passage of energy through the Ryukyu 


Islands. 
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I. INTRODUCTION 


De CouBClTIVE OF THE STUDY 
A serious problem for naval activities, ocean PeauUStELry , 
shore protection and fisheries in the Korean south coastal 
area is the lack of adequate estimation of the waves re- 
sulting from typhoons. The objective of this study is to 
make a partial test of the TYWAVES for forecasting swell 
from tropical cyclones which arrive at a single observa- 
tion site near Cheju-do, Korea (emomiel26. 6-5). 
Forecasting based on wave fields predicted in a 
typhoon area by the TYWAVES [13] developed by NEPRF 
(Naval Environmental Prediction Research Facility) are to 


be verified against swell observations made near Cheju-do. 


pee lYeHOONS IN THE WESTERN NORTH PACIFIC 
During 1979, the Western North Pacific experienced 
28 tropical cyclones. Table I, from "1979 Annual Typhoon 
Pepene 2), shows the significant tropical cyclones for 
that year. Table II from [2] shows the monthly distribu- 
EHOMme Gt tropical cyclones in 1979 and other statistics. 
Most typhoons occur in the summer season and their 
tracks can be classified as one of three typical tracks. 
The first type is that passing south of Taiwan toward the 


west, the second is that crossing over Korea through the 


al 





RyukyuIslands, and the third is that passing east of 
Ryukyu Islands. 


I selected one typhoon of each type; one each in July, 
August, and September in 1979. They are Typhoons Hope, 
Irving and Owen shown in Figures 1, 2, 3 and 4 from [2]. 
tropical cyclone "best track" information is shown in 
Tables III, IV, and V for Typhoon Hope, Irving and Owen, 


respectively. 


C. TYWAVES MODEL (THE TYPHOON WAVES PROGRAM) 

The SOWM (Spectral Ocean Wave Model) run at FNOC 
(Fleet Numerical Oceanography Center) utilizes a coarse 
operating grid system (100~190NM) which does not allow 
sufficient resolution to describe adequately the resultant 
sea state in typhoon areas. Thus, TYWAVES, an improved 
model for typhoons with a locally finer grid, was developed 
by NEPRF. FNOC judged that TYWAVES is worthy of evalua- 
tion aS a possible operational typhoon sea state model 
aay}. 

The TYWAVES is intended to produce fields of signifi- 
cant wave height and spectral wave properties on a mesh 
size consistent with the scales of tropical cyclones and 
is designed, primarily, for the application to the 
western North Pacific. 

The detailed outputs of TYWAVES are in the form of 


Mmolacee at cach of 12 points, of spectral energy 


eZ 





components, Significant wave heights, maximum wave periods, 
and the predominant wave directions, for 00, 12, 24, 48 and 
72 hours, where 00 hours is the time of the first typhoon 
warning issued. An example of those outputs for Typhoon 
Owen 1s shown in Appendix A. 

Using these spectral energy components at selected 
source points, I was able to make predictions based on the 
propagation of these components as swell into the region, 
south of Cheju-do, Korea (see Map I, page 21). 


More details of TYWAVES are described in Ref. ll. 


D. FORECASTING AND VERIFICATION 

Extremely high sea states are known to be generated in 
the quadrant to the right of the direction of movement of 
typhoons. The wave generation in that quadrant of the 
typhoon derived in TYWAVES from a spectral model utilizing 
the Pierson-Moskowitz (1964) spectrum. The model describes 
the spectroangular components of the waves present at a 
Mumoer of grid points in the region of strong winds. Each 
spectral component of interest is then permitted to 
propagate at its appropriate group velocity to the fore- 
cast site. The method is applied to three western North 
Pacmenc EYPneens in 1979 and the forecast products are 
compared with the observed swell data from the south 


coastal area of Cheju-do, Korea. 
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TAGEE |. WESTERN NORTH PACIFIC 


— Sam, 


1979 SIGNIFICANT TROPICAL CYCLONES 


en 


naan CALENOAR MAX WIN NUMBER 

; ; OAYS OF SEG 08S OF DISTANCE 

CYCLONE TYPE NAME OF WARNING WARNING WINO SUP WARNINGS TRAVELLEO 
01 TY ALICE Ol JAN-14 JAN 14 110 930 51 2297 
02 TY BESS 20 MAR-25 MAR 6 90 958 21 1804 
03 My CECIL 1] APR-20 APR 10 80 965 40 2535 
04 i O0T 10 MAY-16 MAY 7 40 984 24 2876 
05 TO TO-05 23 MAY-24 MAY é 30 998 6 2170 
06 TY ELLIS O} JUL-06 JUL 6 85 395 Ze 1612 
07 TS FAYE O01 JUL-06 JUL 6 40 998 20 1837 
08 TO TD-08 24 JUL-25 JUL 2 20 1004 5 1264 
OY ST HOPE 27 JUL-03 AUG 10 130 898 33 3928 
10 is GOROON 26 JUL-29 JUL 4 60 980 13 1058 
1] TO TD- 11 03 AUG-06 AUG * a) 997 14 1088 
le TY IRVING 09 AUG-18 AUG 10 90 954 38 cic 
13 ST JUCY 16 AUG-26 AUG im 135 687 39 2502 
14 TO TO-14 18 AUG-20 AUG 3 20 1006 9 605 
15 TS KEN Q) SEP-04 SEP 5 60 985 13 1418 
16 TY LOLA O2 SEP-08 SEP 7 90 950 23 1298 
17 TY MAC 15 SEP-24 SEP 10 70 984 aS 1831 
18 iS NANCY 19 SEP-22 SEP 4 45 396 14 $28 
19 TY UWEN 22 SEP-O1 OCT 10 110 918 37 2151 
20 TS PAMELA 25 SEP-26 SEP 3 45 1002 6 984 
2] TS ROGER 03 OCT-07 OCT 6 4§ 985 16 1920 
22 TY SARAH 04 OCT-15 OCT 4 110 329 43 1194 
23 ST TIP 05 OCT-19 OCT 16 165 870 60 3972 
24 ST VERA O02 NOV-07 NOV 6 140 915 23 1868 
25 TS WAYNE 08 NOV-13 NOV 6 50 990 22 Teer 
26 TO TD-26 Ol DEC-02 OEC 2 30 998 6 1070 
27 TY ABBY O1 DEC-14 OEC 14 110 951 52 4044 
28 TS BEN 21 OEC-23 DEC J 60 990 10 2245 

1979 TOTALS 14y* 695 


“OVERLAPPING OAYS INCLUDED ONLY ONCE IN SUM. 


TABLE 2. 
1979 SIGNIFICANT TROPICAL CYCLONE STATISTICS ; 

3 

WESTERN (1959-78) 3 
NORTH PACIFIC JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV OEC TOTAL AVERAGE 
TROPICAL ; 
OE PRESSIONS CMO cOl 0 Malo O° | VO ROMEO gi | 5 4.8 
3 

TROPICAL STORMS 0 0 0 0 1 0 2 0 4 1 10 10.0 f 
TYPHOONS Ohm wel Oe *20-- eee Yee 13 ig.0 6 f 
ALL CYCLONES oe [3 “Once 4 cee ee CG 28 32.8 
Q 

8 

| (1959-78) AVERAGE 0.6 0.4 0.6 0.9 1.4 2.) 5.2 6.8 6.0 4.8 2.7 1,3 32.8 : 


FORMATION ALERTS 23 of the 27 (85%) Formation Alert Events developed into tropical cyclones. 





5S of the 28 (18%) tropical cyclones did not have a Formation Alert. : 
: WARNINGS Number of warning days: 149 ‘ 
Number of warning days with 2 cyclones: 38 : 
i 
i 


Number of warning days with 3 or more cyclones: 5 | 
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Table 3. Best track SUPER TYPHOON HOPE 
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O?eauo? lU.2 le?iw en v.49 en UV. -f, Uv, 0.0 awn 0, Ue Me avo Qed de 20, 06 9.9 0.0 Su bao. ne | 
Q7ealeg 10,3 146.4 ¢0 0.0 fen U, <0, Uv, Gn hon 80, “Us Ne Aa ao ie seus 0, oes Sisal. as 
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The dispersive arrival at the prediction site of the 
waves propagating from a single source region is repre- 
sented by a plot of period versus arrival time in Figures 
9, 10 and ll. The energy in each low frequency (7~19 sec) 
spectral component of the swell at Cheju-do is calculated 
by modifying the energy spectrum at the source for the 
Beeraction and shoaling. ~The lowest»period (fT = 4 Sec) 
1s assumed totally eliminated by attenuation. 

The total propagated energy at any given time is 

estimated by summing all the directional components. 
The predicted swell heightsH1/10 are assumed to be related 
to the energy by H1/10=KVE, where K is constant (0.4) and 
E is energy ina unit of ee kiloerg/em- aie hy, Ones eam 
meters. 

By graphing the energy associated with the 5 mean 
period bands as a function of arrival time at the forecast 
site and summing, then a plot of the total propagated 
energy is obtained (see Figures 13, 15 and 7). “has 
plot is compared with the observed wave heights in 


Chapter IV-D. 


E. TERMINOLOGY 
To make clear the notation used in this study, the 


following terms and symbols are defined. 
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sea 


Swell - 


H 1/10 


Heelys 


tar 


h (d) 


Those waves found within a generating area. 
Those waves which have moved out of a 
generating area. 

The distance from the source to the "bottom 
line" where the corresponding periods feel 
“bottom” i.e., the swell travel distance 


in deep water. The depth at the "bottom 


Lo 
—. 


The distance from the "bottom line" to the 


line” is h = 


forecast site. 
AMS Isla Jeyovesverel shel ee ejeceieree ll leew, (i Ss) - 
The average heights of the 1/10 highest 
waves. 

The average heights of the 1/3 highest 
waves (Significant wave height). 

The swell arrival time (GMT). 

The swell travel time. 

Wavelength (Lo in deep water). 

Group velocity (Cg. in deep water). 
Phase speed (Co in deep water). 

The total swell travel distance (D = R_ + R). 
Acceleration of gravity 

Wave angular frequency (2Z/L). 
Wave number (27£/L). 


The ater depth. 


Ze 





O - The angle between wave crest and bottom contour 


O 
at the depth, h = Lo . 
Ty =e Me mcOtalwchengy at "grid point 4 Ineenegenerat— 
ing area. 
D, = he dyreecetonalvenengy Component atbequiadspoune 4 
in generating area. 
Wave spectrum (energy spectrum) - The distribution of 


elther wave height or energy with period. The 

potential energy of the sea surface is proportional 

to the mean of the square of the wave height. 
Dispersion ~- A process which leads to longitudinal 

spreading of the wave energy as, in deep water, 

energy in each spectral component is propagated 

at a characteristic group velocity Cg, the long 


waves having the larger group velocity. 


9 - The mean directional bands of 16 unit points 
rose. 

HS - The wave height before refraction in deep 
Water. 
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II. PREDICTION OF TROPICAL STORM WAVES BY TYWAVES 


Poe DOCATION OF THE SELECTED POINT-SOURCES 

The TYWAVES program produces a printout of the complete 
energy spectrum for 12 points around each typhoon. Each 
point lies at a chosen distance from the storm center in 
a fixed direction, no matter what the storm movement direc- 
tion 1s. These points are representative origins of the 
wave energy emerging from the typhoon. 

Only those points are selected which are in the region 
of maximum winds and where important energy components are 
directed toward the distant prediction site. 

The arrangement of points whose wave spectra are used 


in this study is shown below. 


+9 


+5 +1 +6 


+11 +3 +$ +4 +12 
or aoe +8 
+10 


' phe distances from the storm center to points l, 2, 3 
and 4 are three grid-lengths (3 x 40 = 120 NM) .to 9, 10, 


Maral? are five grid lengths (5 x 40 = 200 NM), and to 


oe 





Syeo,) ? and 8 are 3 x /2 grid lengths (3 x /2 x 40=170 NM), 
respectively. 

Among the 12 points, only points 4, 6, 8 and 12 were 
selected as sources, because the spectral energy from only 
these 4 points can propagate to Cheju-do, Korea (eee N 
MAO Er) . 

The locations of the sources for each time of interest 
for the three storms are shown below in Table VI, which 
was derived and calculated from the postanalysis "best 
track" data shown in Tables III, IV, and V [2]. The 
relationship between latitude and longitude is: 


Lat. = Long. x cos (Lat.) 


iI 


Long. x cos{(storm center + forecast)/2} 
aie 
e.g. : 
When the storm center (best track) is 20.0°N cet ee 


O 


point 4 is same lat (20.0°N), but 120 NM (2~- in lat) 
east of the center. 
nis, 2° Lat = °Long x cos (20-0 ee) 

= Long x cos 10° 

oe Long = De Lat/cos MO = OO or 


26 ts ctee elo 


128.8° E 


iI 


Therefore, point 4 ~* 20.0°N iene =. 
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TABLE 6 


Typhoon Centers and Selected Points Location 


pra (N-PF) 3 =pre (N-PE)) =o pre (CN-°E) = P12. (CN-PE) 


Center (Cn-PE) 


Name Date (Gv?) 


L7s4=l 35-8 
16. 6="3220 
ORG =a 
20.6-126.8 
Zi Sale, 
20.0-130.2 
2220-126 
2350-1236 
24.6-128.2 
2529-128 -0 
2725-1 > 
oo. Odie 7 ae 
31. 7-127 36 
21. 3-193.2 
226-1 
23.5-13273 
24.4-133-3 
Zoomer. 4 
26. oo 
27. 3S. 
2). =e 
206 Doe. 
29.8-134.4 


1534-13329 
16, Ga 3 
Ni eG 28e0 
GsG— ors 
19s 7422 42 
1S. Os1265c 
20-0126. 
Zo L232 
2276-12 6n) 
236 Jat Zo. 
25), Gal 2a 
27. G=1.26.0 
291-1260 
19.8=-132.4 
207 Oris ls? 
Zleo= iol 4 
2224-1330 
Zao 
243 5=132, 0 
253 boca 
25.6-15220 
26. 5-324 
2070 -1oe 


To v4 =133,9 
PAVESI SIS) 
ZIG Z SS 
22 56-L255.4 
2505-0225 
Zan0- 26 nc 
24.0-123 52 
25.5527 0c 
26.6— Zon? 
27. 9-2GeS 
ZOD 26n0 
oo Zon) 
Bea aon 
23a. 
2476-1 o lw 
2555-13154 
26.4-131.6 
Zino woo 
2 on oma 
Coes ec 
ZO oS 2 
BOro=132.4 
Bom seo 


ed 3362 
es Gib oile 
oe Gmileaas 
END (oye Psa 5) 
2s o-l Zee 
20.0-128.8 
OA 
oe ad 
2450-26087 
25,9-126,5 
Zico eon) 
29 6-262 0 
Bll 2G 
Zi 3= 13264 
Zee Oley, 
23. 5-134 
24.4-131.6 
2565-1 te 
26. 5-\3re 
27 6 oo ie 
i Bc: 
25. 5o-L3ces 
Zone lon. 9 


lh) 4S13 8 
TeeG—29.4 
e262 
Z2056-123.2 
Zleoq ZO. 1 
ZOeO=126. 7 
2270-26. 0 
Poe 25.0 
24.6-124.5 
2595 124,3 
2025-12357 
276-123. / 
Seo kee 
Zo 130.3 
PaO 29.5 
Zona Eo e2 
24.4-129.4 
eo) 
2605-02928 
Ziel 2956 
eo b29). 
Zoe om s0... 
2976-13020 


073012 
073100 
073112 


HOPE 


080100 
080112 


IRVING 081300 


081312 


081400 


081412 


081500 
081512 


2 


081600 


081612 
092506 


OWEN 


092518 
092606 
092618 
092706 
092718 
092806 
092818 


092906 
092918 





B. SEA-SPECTRA AT THE POINT SOURCES 

The description of the wave fields using the spectral 
component method developed by NEPRF was used in this study. 
Only those directional components able to propagate to the 
forecast site need be considered. 

A schematic example of the TYWAVES spectral prediction 
in a typhoon area is given below. The whole spectral table 


for the selected sources is shown in Appendix B. 


T - 8 spectrum 
(80 components) 





ia omee cal 


"5 ~ mean | 
‘period bands 


Table VII is one example of a period-direction spectrum 





at source point 4 in Typhoon Owen, at 06 on 25 September 
S79. 
The mean period is 13 sec with the maximum energy 86 


(2° 


. kiloeng/em*) in the dominant period. The dominant 
wave direction is SE, but the energy from that direction 
will not arrive at the forecast site, only the SSE component 


wiil arrive at the site (as seen in Figure 5). 
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C. SELECTION OF THE DIRECTIONAL ENERGY AT POINT SOURCES 

A constant energy was assumed between parallel rays 
separated by width 80~160 NM (2-4 grid length), which is 
40~80 NM right at the grid point depending upon the typhoon 
size. This is why the energy in each component is assumed 
constant along the deep water propagation path. As example 
of possible refraction diagrams under the above considera- 
tion and bottom contours in Map 1 is shown in Figure 5. 

Only the SSE component was considered to propagate to 
the forecast site. The SE component would refract toward 
China, west of Cheju-do, and the S component would pass to 
the east of Cheju-do as seen in Figure 5. 

All the directional components were chosen in the same 


way, and are underlined as shown in Table 7. 


Figure 5. Possible 

nN refraction diagram at 
point 4 at 06 GMT 
September 25, 1979. 
The curves around 


aes 
A i ‘| 
\ od 7 
x | EVONoon center are 
fi Ss | contours of wave height 


Oe (H 1/3 in feet). The 
\ 


aetauled bottom contours 
1OUFM «ci | 
LINE 





are shown on Map l. 
aa: 


+ . Tr = - 
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TABLE 7 


Period-Directional Spectrum at Point 4 at 06 GMT eepeeu ee” 
25, 1979. The energy units are 2°- kiloerg/em2. 


TYPHOON NAYES 20 ee 


— ee ee se a ee ee 


Haus = a5 14 


NORTH-CORRESPONDS TO THe AZ ie see 


N OS U i J : : i 

NNE J) ge ¥ S a : : Is 

ae Je e G : s a C 

ae cs , - s : 4 

= he u U B ; g J 

as Se J 6 ers 1 U 

So Gout o 6. eae 2 3 

r Sie alae J Ces oes 3 L 
> es L a | ls 2 

SSA 0 U : 2 2 C 

= SA ag U s : L g J 
WSW G85 C v L : C c 

W one J g g E u J 

aa “WN ora 9 ¢ 3 S C c 
NW Bae ¥ C v : 5 C 

NNW a . L ts . C id 
a a Se ee 

TOTAL 232 


J 
PERIOOS 4 


Meeetousn AND DIRECTIONAL ENERGIES AT POINT SOURCES 

Figures 6, 7 and 8 show the energy components at each 
selected grid point in the typhoon area. Each line is 
labeled by selected grid point (Source point). ae indicates 
total energy given by TYWAVES at point n and Dd, represents 
the sum over all periods of the energies with proper 


feeeecer1on tO reach Cheju-do.. 


a0 





D(8.) —— E, (T,, 9.) 
For example, where n = 12, 8 is SE, then 
D)5(8=SE) = EL (T=4, ao) eee LO), 6 0 =S8E) 
pee, (T=13, Oto) Ee epee sr) 


~ E,9 (T=19, OSE). 


This represents the total energy at point 12 directed 
from SE toward Cheju-do. 
The energy 1S in practice converted into height H 1/10 
by making H 1/10 = 0.4 VE, which is described in Chapter I.D. 
some special features of the wave fields in each source are 
now discussed. 
1. Typhoon Hope 
The maximum sustained wind speed was over 100 kts 
for two days from O0OGMT July 31. The maximum wave height (H 1/10) 
recorded was 12.80 m. The directional energy components from 
grid point 12 were ignored from 00 GMT July 31 to 00 GMT 
August Ol, because they were too small. 
Paes yonOoOn Irving 
Only the propagation of energy components until 
00 GMT August 16 was considered, because the forecast site 
1s inside of typhoon area after that time. Therefore, after 
00 GMT August 16, the propagated wave energy should be 


combined with local energy at prediction site. 


Si 





Among the directional components, only those from 
grid points 4 and 12 were propagated to the prediction site 
during the typhoon period, because the energies from other 
sources were refracted away from Cheju-do. 

3. Typhoon Owen 

After 06 GMT September 29, the typhoon passed east 
of southern Japan. No further propagation to Cheju-do 
Seecunrred. 

In most cases, only one directional energy component was 
selected at each grid point. But both components (SE and 
SSE) were able to reach the forecast site on 18 GMT September 
27, from grid point 4, because of the width of their ray 


Boundaries at distances near the source. 
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Peek OreLrON OF SWELE CHARACTERESTICS 


The periods of high waves due to typhoons are usually 
in the range of 7~19 sec. They travel in deep water with 
the group corresponding velocities 10~30 kts. It normally 
necessitates 1~4 days for conspicuous swells generated in 
a typhoon in the far south seas to reach Cheju-do, Korea. 

The group velocities of periods 7 and 10 sec are 
constant along their paths since they remain in "deep 
Poets 2G for periods 13, 16 and 19 sec, the group 
velocities vary because of depth variations. Thus, I 
considered two distance components: in deep water R, and 
in shallow water R for periods 13, 16 and 19 sec. 

The travel distances between grid points and observa- 
tion points are calculated with the assumption of plane 
geometry. Then, group velocities, travel distances and 
travel times are derived for each period in this section. 
Peed SEETENG OF THE FORECAST SITE AND ITS RELATION TO 

THE WAVE SOURCES 

The passes of the Ryukyu Islands act as windows between 
the energy sources and the forecast site, Cheju-do, Korea. 
Since the largest land length for blockage of energy is 
enayeabeut 40 NM, this length is not Sufficient to inter- 
mipe totally the energy propagation. Thus, as a further 


simplification, I neglected the effect of the Ryukyu Islands 
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against the energy propagation. The window is assumed 


aeemeewent FOr total propagation. 


B- GROUP VELOCITY 

Some Simplications are used in assigning group velocity 
values along the route from typhoons to the observation 
Site. 


From linear wave theory group velocity, e 1S given by 


2kh 


a! __2kh 
Same |S * Sten ob. 


= ee 


(all the symbols are defined below) 


Z tanh kh. In deep water 


and phase speed c is given by c 
this is well approximated by eee G/G = 32050  “Cebe): for 

T given in seconds, and in shallow by c = Ygh. But in the 
general case ¢ varies with both the depth of water and the 
wavelength. 

The classifications "deep" and "Shallow" are given in 
terms of "relative depth h/Los described below. Since 
there is no "Shallow water" between the wave sources and 
Cheju-do, Korea, (corresponding to the periods used by 
TYWAVES), I will use only the general equation and the 


deep water approximation. In these equations 


Cy = group velocity 
c = phase speed (co in deep water) 
nm = the fractionof energy propagated at phase 


speed 
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g = gravity acceleration 


h = water depth 


G6 = wave angular frequency (2T/T = 2kf, where T is 
period) 

k = local wave number (2T7/L, where L is wave length) 

£f = wave frequency (1/T) 


Poromimportant to note that c and group velocity, Cg 
Moseeoe LOund aS functions of both h and o. A straight- 
forward method is illustrated below, where the dependence on 
kh is replaced by the deep water wavelength Ly = te c/a. 

The local wavelength is then found from L/L: a function 
of the relative depth h/L: 


Since C/C. = L/Los the phase speed c can be found from 


head O . 
To estimate the group velocity by definition, Ca = 
1 2kh os ; 
5 C (1 + eaaaorn = nS where, n takes on the following 
values for corresponding h/L,} 
1 1 1 
general case (= > h/L, Ei 59)? lon>> ; 
deep water case (h/L Se n = £ 
Op = Dev 
In order to simplify calculations, I used n = 3/4 in 


the waters between the Ryukyu Islands and Cheju-do 
where, for the longer periods, the relative depth range 


5 >h/L>s5 . With this simplication, the corresponding 


group velocities become the following: 
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T = 13 sec; oe = 29 45 ts 
T= 16 sec; Cy = 33.6 kts 
T = 19 sec; oe = 35.8 kts 


The influence on the calculation of ae of the use of 
these representative constant values is discussed in 
@iapter IV.D. 


For deep water 


iS 


Bs SUS (kts: for £ in sec) 


1 
Soe” nc = 5 Ce =), Seo ees wom — sn See) 


Therefore, group velicity in knots and depth at which 


C_ replaces C are shown below: 
g So 


T (sec) 7 10 1S 16 i 


a (Co iiss meme Oi ge a2 Cee 2G) nares acy 2 Od 
O 


L 
h(+- in ft) 126 256 433 656 924 


C. TRAVEL DISTANCE AND ARRIVAL TIME 
The travel distance between grid points in the storm and 
the observation point 33.2°N 126.6°E are calculated with the 
assumptions of plane geometry. Consider this example. 
On 12 GMT July 30, the grid point 4 was at 17.4°N133.90R8. 


Thus, the north-south component is given by 
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K 
il 


lecmrsute)> = lat Corea point) 7] x 60 NM 


il 


[eee2 = 17.4; x60 


948 NM 


and the east-west component is approximated by 


i [long (grid point) - long (site)] cos (mean lat) ] 


x 60 NM 


peer 126-6] 6 cos (17-4) eGo 


418 NM 


Thus, the travel distance, D, is given by 


1. 
DEEMC ety-)\* and for guia poune 4, 
1. 
D = (418% + 948°)% = 1036 NM 


The travel time, t, is given by 


ee R aa 
t= VY + E , where Ra + R= D 


Ro = the deep water distance 


R= the remaining distance 
iiiererore, the arrival time, tar of each period band is 
given by 
eve (Re £, + t, where t. = leaving time from 


typhoon area 


i travel time 
The Figures 9, 10 and 11 show the swell arrival time 
hOr C€ach of the three typhoons. Each curve indicates the 


arrival time of energy from a source at a specified leaving 


time labeled with various periods (from 7 to 19 sec.). 
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Appendix C shows the arrival time calculations for 


each typhoon. 
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IV. THE ENERGY PROPAGATION PROCEDURE 


A. PROPAGATION OF SPECTRAL ENERGY COMPONENTS 

The spectral forecast permits the forecaster to deal 
with only that range or periods which have important energy. 
Each spectral component is tracked with its respective group 
velocity, only those directions being chosen for which waves 
can reach Cheju-do. When the typhoon is moving with a 
velocity component toward Cheju-do, only spectral components 
with a group velocity greater than the movement of the 
typhoon are considered (as in Typhoon Irving). 

Wave energy generated in a relatively small region at 
all frequencies will spread over a much larger region as it 
propagates outward from its source, and the wave character- 
1lstics change in such a way as to become more "Swell like." 
There are essentially three processes which contribute to 
this change in wave characteristics: dispersion and angular 
Spreading, which are modified by shoaling and refraction, 
and attenuation. 

In the procedure followed here, I have ignored attenua- 
tion for waves of long period. Dispersion and angular 
spreading are accounted for by simply following components 
to Cheju-do with appropriate shoaling and refraction 


factors. 
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Por ease of calculation Shoaling and refraction processes 
are simplified as described below. As seen on Map l, a 
Slightly curved contour of 100 fathoms (590 ft) connects 
the northern tip of Taiwan to the south-western Elp Ce 
Japan. The slope from 590 iG (ice Inas io ells Si (70 
fathoms) is very steep. But most bottom topography along 
the path to Cheju-do;, north of the 413 ft bottom contour 
is almost flat with the depth of 45 fathoms (266 Pe) pmo 
the forecast site (33.2°N 126.6°R). 

Therefore, I have approximated to underwater topography 
by assuming only two water depths, a deep water and a 
shallow water (intermediate water depth) region with an 
abrupt jump between them. 

To compute the wave characteristics at a shoal water 
Bice, the Shoaling and refraction are considered using 
the values of oe and n from Chapter III.B. 

The energy of component EY (T,9) in the typhoon area 
(deep water) is transformed after shoaling and refraction 


to its value at the forecast site ACCOEGING “EO 


2 
e 


E(i, j) = EL (i, j)-k* G4, 9) k? GG, 4) 
O Ss 

Wieme eB o(i, 7) is the energy of the component of 

period le and which had the direction 8, in the 

generating area, 

and Ne (te. 3) 5 Ky (1, j) are the respective 


Shoaling and refraction factors of those components. 
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Petatis of calculation of Ke and K are given in the 
following sections. The refraction and Shoaling calcula- 
tions for three typhoons are shown in Appendix D. The 
following is a sample calculation on 12 cm? July” 3UMeen 


Typhoon Hope at source point 4. 


it 7 10 13 16 19 

K 1 1 0.92 0.89 0.90 

kK. 1 1 i 0.989 02972 
(Ke ee 1 1 0.8464 0.7748 0.7653 

E, 7 19 wlth 9 14 

EB 7 19 9.3 7.0 neha 


where Ke isederived Grom plate G—l(6] 
with h = 40 fathoms = 236 ft 
KS is derived from Figure 2-19 [5] 
with h = 413 ft and a, = 40°. 


B. SHOALING AND REFRACTION FACTOR 
ieeoloaling Factor 
The shoaling factor, Ke = H/F} was derived from the 
plate C-1l [6] with depth and period of the appropriate 
Spectral component. The shoal site depth near Cheju-do 
(33.2°N 126.6°E) is 40 fathoms, and the relative depth h/L, 


and shoaling factor of each period are shown below. 
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pl ae 10 diigs 16 I, 


A/L, 0.94 0.46 Oza 0.18 0.13 
ie i il: 0.92 0.89 0.90 
Po KecEracuion Factor 


Generally, two basic techniques of refraction 
analysis are available--graphical and numerical. Several 
graphical procedures are available, but all methods of 
refraction analyses are based upon Snell's law. Refraction 
teayeee treated analytically in any region with straight 
and parallel contours, by using Snell's law directly: 


; Cc . ; 
sina= a) sin a), where a is the angle between the wave 
O 


crest and the bottom contour, and a, is the angle between 


O 


the deep water wave crest and the bottom contour. 

Figure 2-19 [5] shows the relationships among 
or Woe period, depth and refraction factor in graphical 
form. I derived the refraction factor from using this 
Grapm, Figure 2-19, the bottom contour and period. Also, 
I assumed that the refraction factor is 1 if the angle Cl 
(between crest and bottom contour) is less than 10° and 
that the refraction occurs only one time at the depth of 
maximum bottom gradient where waves are refracted because 


of the wide flat-bottomed portions over most of the inter- 


mediate water propagation path as seen on Map l. 


50 





With these considerations, I derived the refraction 
factor and predicted the shoaling energy. Those procedures 


are shown in Appendix D. 


C. SHOALING AND REFRACTION OF THE SPECTRAL COMPONENTS 

Each deep water wave spectral component derived from 
the TYWAVES model moves with its respective group velocity. 
I considered the 80 energy components to behave as mono- 
chromatic component waves. To assess the shoaling and 
refraction of each to the shoal-water site I assumed that 
the wave power transmitted between a given pair of 
orthogonals remains constant at all depths (this means no 
frictional losses, diffraction or scattering, and also 
implies that a steady state exists). 

With these assumptions the wave power P is given by 


Pee Db = EF. C b 


g ee 
oe bo 
_ O 3 ee es 2 
Thus, E = EG oe ES Be Ks Ke = oe PgH 
g 
where E = the average wave energy per unit area 


of sea surface for waves transformed by shoaling and 


metract1On). 
P = water density 
g = acceleration of gravity 
H = wave height of transformed waves 
GC Cue group Vvelocit 
G5’ 9 g Pp iy 


Sil 





or 
oy 
il 


orthogonal separations 


Ww 
Il 


shoaling factor 


nw 
I 


refraction factor 


The energy in each spectral component of the swell at 
the observation site was calculated by modifying the energy 


spectrum at the selected sources for the effects of shoaling 


Ziiyiee 
S a 


and refraction according to equation E = E,°k 

The energy associated with the various components as a 
funct1On Of time of arrival at observation site, as seen in 
Begures 9, 10 and 11. 

The total energy in the swell at any given arrival time 
is estimated by summing all the shoaling components at that 
time. 

In summary, each component in typhoon area's T - 6 
spectrum is shoaled and refracted using the i and kL 
values appropriate to it to find the energy at the observa- 
mon Site. 

1. ShoalingProcess (Computation) 

See Appendix D. 


2. Shoaling Energy Components From Each Source Versus 
Arrival Time 


As seen in Figures 12, 14 and 16, the energy of 
the components from all sources is shown as a function of 


its arrival time at Cheju. See Appendix D. 
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3. The Predicted Swell Waves at Prediction Site 

As shown in Figures 13, 15 and 17, the predicted 
swell waves are the sum of all transformed components at 
given time at the forecast site. 

For Typhoon Irving (see Figure 15) beginning 
12 GMT August 16, the forecast site is already inside 
typhoon area. Therefore, there is no prediction done 
aeeer that time. These predictions are discussed in 


Giapter IV.D. 


De THE OBSERVED DATA 

Observations of wave conditions for verification of 
the swell forecasts were obtained from the sources, the 
National Oceanographic Data Center (NODC) and the Republic 
of Korea (ROK) Navy. All listed heights are based on 
visual estimates. 

Following Table VIII shows the visual observation of 
swell heights (H 1/10) for each typhoon. The data from 
NODC are sparse and often far from the forecast site, 
but they appear to be samples from the same set as those 
of the ROK Navy. The data from ROK Navy visual observa- 
tions were made at Be ake 126.5°E close to ty Od Memo r 
interest ere N, 126.6°E). So I considered this Somme 


exactly the same as my forecast site. 
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TABLE 8 


The Visual Observation Data by NODC and ROK Navy 


Pee Or bypnoon Hope 
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DTG (GMT) Location GN ©) H 1/10 (meters) Source 
79080100 Dor ~ ee Saye 25 NODC 
0104 B25 — 26.0 2e5 NODC 
Galele2 Bie 2-167 4.0 NODC 
Orla3 Zilicr) > Ze..6 Srno NODC 
121 Cor 25. Gigs NODC 
0200 Zoe = 126.0 6-5 NODC 
0206 Ze= 0-129. 0 Se NODC 
0218 Gr orl 5G 520 NODC 
0300 28. 1-12 8..6 Si NODC 
0306 31.4-126.2 320 NODC 
Oe12 S25 125.9 Eas. NODC 
080100 Boe2— 26. 5 eau ROK Navy 
0104 33.2-1265.5 eo ROK Navy 
0108 BGe2—1256.5 On ROK Navy 
Ourie2 Be. 2-126. 9 4.0 ROK Navy 
0116 Ben2- 26.5 4.5 ROK Navy 
0120 Bees = loo. 5 4.5 ROK Navy 
0208 Bon 2-126. 5 5.0 ROK Navy 
0212 Bie 2 26> a5 ROK Navy 
0216 See 126.25 55 ROK Navy 
0220 eee 126.5 a5 ROK Navy 
0308 8362-126, 5 Ba ROK Navy 
em2 Borne 12665 Ge ROK Navy 
Geis See Zoe 520 ROK Navy 
0320 33.2-126.5 Bec0 ROK Navy 
0408 83.2=126..5 4.0 ROK Navy 
0412 BSa2— 2655 Be ROK Navy 
0416 3322-26. 5 520 ROK Navy 
0420 Some = leo Dees ROK Navy 





For Typhoon Irving 
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DTG (GMT) Location (282m) H 1/10 (meters) Source 

79081400 29 O12 6 a NODC 
1403 28.5-127.6 4.0 NODC 
1406 28.7-129.3 345 NODC 
1409 Cie = eo] 4.0 NODC 
1412 Doe 2a loads 1 4.5 NODC 
S00 B00 —12 5. 0 G0 NODC 
135.03 B07 Shen 3 4.0 NODC 
les) ils Zoo oU. 3 G0 NODC 
ieee 24 2 1 30%.0 9.0 NODC 
iLSy-ANl 2655-26. 0 SU NODC 
1G0'0 27 = 130 2.0 Grae NODC 
1603 Dio ao.) 8.0 NODC 
1606 Pee 12 9 2 Seo NODC 
i700 eon lo <5 Gn NODC 

79081400 Bone 126. 5 oe ROK Navy 
1408 bee 2-120. 5 Des ROK Navy 
1412 Beel 126. 5 S40 ROK Navy 
1416 Riga = 120. 5 So ROK Navy 
1420 So 2-126% 5 Psd ROK Navy 
1508 Bee -i 26.5 4.0 ROK Navy 
sal Bee = 6262.5 4.5 ROK Navy 
ILS ks: Boece 20.2 Sue ROK Navy 
iS 2.0 Bene 2645 6.0 ROK Navy 
1608 Biiga= 126.5 Sit0) ROK Navy 
oi 2 Be 2 = 2 65 9.0 ROK Navy 
616 Biome = b26. 5 ne ROK Navy 
1620 835.2 = 126.5 8.0 ROK Navy 
1708 8362-126, 9 oa0 ROK Navy 
hy eZ 2382527 020.5 9.0 ROK Navy 
LG BOmee b26. 8.0 ROK Navy 
20 Some al26. 5 a0 ROK Navy 
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Er COMPARISON OF SWELL PREDICTIONS AND OBSERVATIONS 

The predicted and observed heights are plotted in 
Figures 13, 15 and 17, and are good agreement. Detailed 
comparisons for each typhoon follow: 

1. Typhoon Hope 

The first predicted swell arrival times from each 
of three sources, and the NODC and ROK Navy observations 
are almost the same with similar energies (see Figure 13). 

The peak energy arrival times are also nearly the 
same, but the peak predicted energy is slightly higher than 
the observations. This is reasonable agreement since it is 
estimated that the visual observations can be error. Verploegh 

1961 [18] estimated the average observational error for a 
visual observation of wave height varies from 1 ft at 5 ft 
wave heights to 3 ft at 18 ft wave heights. 

Zee ye neon irving 

As seen in Figure 15, the predicted swell arrival time 
nearly matches the time shown by both sets of observations. 
However, the predicted energy of the rise is slightly higher 
than the observational values. The peak energy (around the 
time of 20 GMT August 15) is also the same as the observa- 
tional peak. 

After 00 GMT August 16, the prediction’site is within 
the typhoon's wind circulation area. Most wave heights (H 1/10) 


are 8m and the highest observation is 10m. These waves 
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dominate the entire wave field and no swell can be dis- 
tinguished later than 00 GMT August 16. 
3. Typhoon Owen 

In Typhoon Owen (see Figure 17), the swell arrival 
time indicated by observations is earlier than the predicted 
time by about 6~12 hours. The predicted peak energy lies 
close to the observational values. 

In fact, the NODC values are higher than those of 
the prediction, and ROK Navy observations are lower. 
PaeGepesmear the peak both observation sets show similar 
time variations. Therefore, the observed values are con- 


Sidered consistent and probably accurate. 


Beak Rkorn, SOURCES 

I have made several assumptions in this study in order 
to simplify calculations. The most serious error sources 
involve assumptions about the windows in the Ryukyu Islands, 
group velocity in shallow water, and simple bottom contours 
in shallow water. 

There are also differences between the predictions (which 
include waves dependent on the local weather conditions) that 
make it difficult to evaluate the prediction. Regarding the 
assumptions about group velocity in Chapter III.B, I used 
n = 3/4 in shallow water for calculation of group velocity 
C =nce. In Typhoon Hope, which has the longest shallow 


g 
wWieeamtrave! distance, miscalculation of or would have its 
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greatest effect. Yet the observed and predicted energy 
peaks are not greatly separated. This may mean that the 
approximations are realistic. 

The bottom contour assumptions seem to be reasonable 
through three tests. The local weather condition is the 
most serious factor contributing to differences between 
the observations and the predictions. The local wind 
pattern was the following at Cheju-do, according to the 
ROK Navy data set. 


1979 072900-073100: SW 5-8 kts (preceding Hope 
swell) 


080100-080300: SW 6-10 kts (during Hope 
swell) 


1979 081200-081300: S 6-8 kts (preceding Irving 
swell) 


Cole GO-031700. 4) 5S 30=60) sts (dieing fay onc 
swell) 


092400-092600: E 5-8 kts (preceding Owen 
swell) 


092600-093000: E 20-30 kts (during Owen 
swell) 


Therefore, before the swell arrivals, the local wave heights 
were conSidered less than 1~1.5 m. During Irving and Owen 
there were important local sea contributions to observed 
height. 

Lastly, the original error sources, i.e., those involving 
input parameters for the TYWAVES model, the location of the 


typhoon center, the winds in the source region, and typhoon 
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size, etc., are ignored in this study. Those errors are 


@eseussed in Refs. [2] and [17]. 
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Vi "GONGEUSTONS 

Pronettemmes 13, 15 and 17, I can conclude that most 
predicted swell heights are lower than those of the obser- 
vations (combined sea and swell). 

The times of occurrence of the predicted peak heights 
agreed reasonably well with those of observations for the 
swell from each typhoon. 

These results on the basis of this limited test suggest 
that TYWAVES predict satisfactorily those swells in the East 
eitaa coca whieh originate in tropical cyclones in the 
western North Pacific. 

Computer aided predictions may improve the quality 
of the forecasts by reducing the need for simplifying 
approximations, as in the case of the treatment in this 
thesis of the shoaling and refraction processes. Such 
predictions would also provide a larger base for assessing 
the accuracy of the method. Additional verification is 


required to draw more specific conclusions. 
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APPENDIX A 


Outputs of TYWAVES For Typhoon Owen 
at 12 GMT September 26 


The following tables and figures show the computer anal- 


ysis for Typhoon Owen by TYWAVES. 


ae 


b= 


Peptod=cdteect1onal spectrum at each of 12 sources. 


The Significant wave height (H 1/3 in feet) distri- 
bution in typhoon area. "O" indicates the land 

area and the distance between the grid points is 

40 NM. 

The maximum wave period distribution in typhoon 
area. "-1" indicates the land area. 

The maximum wind wave directions distribution. 

The direction indicated with 16 unit point rose, from 
Nomriwiel = to uNNE (6) with CeWedirection. 

NEPRF Typhoon Wave Program Analysis for Typhoon Owen 
during the typhoon period from the time of the first 
typnoon warning to 72 hours later, based on post- 


analysis data. 
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APPENDIX B 


sea Spectra at Source Region for Three Typhoons and Significant 


Wave Distribution 


The following tables and figures show the period-direction 


Spectrum at selected source grid point and NEPRF typhoon wave 


program analysis for wave height (H 1/3 in feet) distribution 


around the typhoon center, respectively. 


NEPRF Typhoon Wave Model (period-direction spectrum 
at each selected grid point) for Typhoon Hope. 
NEPRF Typhoon Wave Program Analysis for Typhoon Hope. 
NEPRF Typhoon Wave Model (period-direction spectrum 
at each selected source grid point) for Typhoon 

Ie Valiy 

NEPRF Typhoon Wave Program Analysis for Typhoon 
einvelnia Gh 

NEPRF Typhoon Wave Model (period-direction spectrum 
at each selected source grid point) for Typhoon 
Owen. 

NEPRF Typhoon Wave Program Analysis for Typhoon 


Owen. 
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APPENDIX C 


Propagation Energy Arrival Time 


The following tables show the components from each 
source grid point, travel time and arrival time (GMT). 
Ne and R are determined from Map 1 (U.S. Navigation Chart 


No. 94027, Scale 1: 927.700 at lat 32°15') from the critical 


depth of each period. 
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APPENDIX D 


Shoaling of the Spectral Energy Components 


For refraction factor computation at 00 GMT August 13 
Paoierolit ston Vyonoonwmirvying, the critical water depth is 
325 feet and the angle between the swell crest and the 
Poe eOn cOnmeOUuEemms25 £E depth 1s 40°. Thus, with h/gT* 
and the Figure 2-19 [7], the refraction factors are derived 
like the following: 


T (sec) = 7 EO les 16 19 


=). are L 0293 05968 





In the same way, the other refraction factors are shown in 
shoaling computation, in the following tables: 
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